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Abstract 

A systematic procedure for deriving the master equation of a dissipative system is reported in 
the framework of stochastic description. For the Caldeira-Leggett model of the harmonic-oscillator 
bath, a detailed and elementary derivation of the bath-induced stochastic field is presented. The 
dynamics of the system is thereby fully described by a stochastic differential equation and the 
desired master equation would be acquired with statistical averaging. It is shown that the existence 
of a closed-form master equation depends on the specificity of the system as well as the feature of 
the dissipation characterized by the spectral density function. For a dissipative harmonic oscillator 
it is observed that the correlation between the stochastic field due to the bath and the system 
can be decoupled and the master equation naturally comes out. Such an equation possesses the 
Lindblad form in which time dependent coefficients are determined by a set of integral equations. 
It is proved that the obtained master equation is equivalent to the well-known Hu-Paz-Zhang 
equation based on the path integral technique. The procedure is also used to obtain the master 
equation of a dissipative harmonic oscillator in time-dependent fields. 
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I. INTRODUCTION 



In the real world physical systems are not exempt from the disturbance of their surround- 
ing environment. Because of dissipation caused by system-environment couplings [liJsl, the 
dynamics of open systems will change dramatically from reversible to irreversible behav- 
ior. The most representative model of the dissipative dynamics is the Brownian motion, 
a subject extensively investigated in many academic areas ranging from natural sciences, 
engineering to social sciences U0. As one observes that the environment makes the dy- 
namics of the system a stochastic process, it is naturally to use random noises to imitate 
the influence of the environment, which results in the phenomenological Langevin equation 



and Fokker- Planck equation methods 



The Langevin equation describes the motion of 



a classical Hamiltonian system subjected additionally to a noise due to the environment, 
while the Fokker-Planck equation depicts the evolution of the phase-space probability for 
an ensemble of the system. The most successful microscopic model for the environment 



12| . which has been shown to 



consists of infinite number of effective harmonic oscillators 

be a generic heat bath [isl [l^ . The harmonic-oscillator bath with linear couplings to the 
system under study is called the Caldeira-Leggett model in the literature. 

Traditionally dissipative dynamics is an essential subject of nonequilibrium processes 
and plays an ^jndispensable role in condensed-phase dynamics, transport as well as spec- 
troscopy 



15 



171]. Nowadays the rapid development of quantum information and quantum 



measurement provides extraordinary paradigms for exploring di ssip ative effect, especially 



18 



2l| . Because it is an 



the interplay of quantum coherence, dephasing and relaxation 
inherent many-body problem, dissipative dynamics is usually very difficult to solve exactly. 
There are, however, still several general frameworks available for investigating the diversified 
features of open systems . For instance, the projection operator technique developed 



by Nakajima 



22| and Zwanzig j23l] has been widely used in fields of spectroscopy and quan- 



tum optics where dissipation in general is weak and thereby perturbation treatment can be 
employed. The infiuence functional method based on the path integral, on the other hand, 
which was first proposed by Feynman and Vernon 2J] and popularized by Caldeira and 
Leggett [isl has been shown a powerful tool for theoretical analysis [l|. This method has 
also been implemented as a numerical technique to simulate dynamics of dissipative two- 
state and three-state systems 25|, |26| . Following Feynman, Stockburger et al. interpreted 



2 



the influence functional in the Caldeira and Leggett model as a random field and put for- 
ward a stochastic formulation for the dissipative dynamics Along this line one of the 
authors extended the stochastic idea to general heat baths and was able to establish that the 
influence of the heat bath on the system can be fully characterized by its induced stochastic 
fields 28(1 . Based on the stochastic description, a deterministic approach comprising hierar- 



chical equations of motion was developed 



29 



32| . This scheme is only appliable to specific 



dissipation where the bath-induced stochastic field is the Ornstein-Ulenbeck type or the like. 
Fortunately it is shown, that the hierarchy method combined with stochastic realization can 
be used to accurately simulate zero-temperature dynamics of dissipative two-state systems 
at strong dissipation [29 1. 

Analytically solvable models are always desired because they provide deep insights into 
the understanding of dissipation and benchmark results for comparison when developing 
approximations. It is unfortunate that in dissipative dynamics only few systems are an- 



alytically solvable. Good examples include models of the harmonic oscillator [33N50| . the 
free particle 5l|], and the parabolic barrier 52l-l55l|. In this paper we will show that the 
stochastic decoupling approach can be invoked to obtain the master equation of the dissipa- 
tive harmonic oscillator with or without driving fields in a straightforward and elementary 
way. Although some results were already reported and the method was briefly outlined 



in previous work 



38 



39|, a full procedure by which one can follow step-by-step is yet to 



supply. This paper will fill this gap. We will give a detailed report on the derivation of 
the master equation for dissipative linear systems from the establishment of the stochastic 
differential equation of the dissipative system to the statistical averaging that leads to the 
desired master equation. 

The paper is organized as follows. In Sec. II we briefly review the stochastic description 
of quantum dissipation. In Sec. Ill we apply the stochastic formulation to a dissipative 
harmonic oscillator and derive its exact master equation. In Sec. IV we prove the equivalence 
between our results and the known results, Hu-Paz-Zhang master equation. In Sec. V we 
extend the scheme to the case of the harmonic oscillator in time-dependent external fields 
and work out the master equation. We conclude the paper in Sec. VI. 
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II. THEORY 



A. Stochastic description: Primer 

To study the dissipative dynamics, we start with a system-plus-bath model defined by 
the Hamiltonian, 

H = Hg + Hh + Hint, (1) 

where Hg is the Hamiltonian of the renormalized system of interest, Hf, the Hamiltonian of 
the bath and Hint the interaction energy between the system and the bath. Without loss of 
generality the couplings can be written as Hint = J2a faQa, where /„ are the operators of the 
system and Qa the operators of the bath. According to quantum mechanics the evolution of 
the total system obeys the Liouville equation, namely, 

^h^=[H,p{t)], (2) 

where p{t) is the density matrix. Keep in mind that we are only interested in the dynamics 
of the system. The reduced density matrix Ps{t) = Tibpit) provides sufficient information 
we require. Because the complexity of dissipative dynamics lies in the coupling between 
the system and the bath, it would be desired to decouple the interaction in such a way 
that the evolution of the bath will no longer be explicitly involved in the evolution of 
the system. Actually, it was shown by one of the authors that the Hubbard-Stratonovich 
transformation can be used to convert the system-bath coupling during the evolution to the 
correlation of stochastic processes for the separated but random system and bath [28] . Doing 
this, the price one has to pay is to introduce auxiliary stochastic fields in subsystems. Later, 
it was recognized that a simpler and natural language way to formulate the dissipative 
dynamics as a stochastic one is the Ito calculus [28 1. 

The Ito calculus is concerned with a Wiener process W{t) = J^dt'p{t'), where p{t) IS a 
Gaussian white noise with zero mean and delta function correlation, i.e., M{p{t)} = and 
M {p(t)p(t')} = 5{t — t'). Roughly speaking, the white noise p{t) can be regraded as a series 
of independent random numbers at time slices to = 0,ti, ...,tAr = t, where one can simply 
assume the time steps At^ = tj — tj^i = At are uniform and At — t- 0. Thus the distribution 
function for any pj = p(tj) is 



, , At 
P \ Uj) = lim \ — e 
^^^^ At^o V 27r 



For a random process F[ii\ or F {/ii, ...,iin) in discrete-time representation, the stochastic 
averaging M means the expectation value of F {/ii, /in), i-e., 

M{F[fj,]}^ / l[[dfijP{f,j)]F{f,i,...,f,N). 

Note that dW is of the order ^/di and {dW)'^ = dt. To apply the Ito calculus we now 
consider the Liouville dynamics Eq. (2) with a disentangled initial state p(0) = Ps{0)pii{0). 
When calculating the time derivative of a composite stochastic process constructed from the 
Wiener process, one should take into account the second order contribution of dW. Let us 
analyze the following stochastic differential equations. 



ihdps 



Hs,ps dt+ ^^^fa,ps dWia + i^^ya,Ps} 



dW. 



2a 



(3) 



and 



ihdpb 



Hb, pb dt + — Y^ [ga, Pb] dW2a + ^ {9a, Pb} dWi^, 



(4) 



where Wia{t) = J* dt' [pia{t') + ip4a{t')] and W2a{t) = /J dt' [p2a{t') + ip3a{t')] are complex 
Wiener processes with independent Gaussian white noises Pka{t) (A; = 1 — 4). Apparently 
there is no direct interaction between the two subsystems and the equations of motion for the 
system and the bath Eqs. (3) and (4) are independent. Using the Ito calculus of complex 
Wiener process, dWadWk = dW*dW^ = and dWadWk = 26akdt, and nonanticipating 
property, that is, the combined stochastic process ps{t)pb{t) is statistically independent 
of dWi{2),u M{ps{t)pb{t)dWi{2),t} = 0, we can readily prove 



iM {M {ps{t)pb{t)}) = ihM {dps{t)pb{t) + Ps{t)dpb{t) + dps{t)dpb{t)} 



Hs + Hb + J2 f'^da, M {ps{t)pb{t)} 



dt. 



(5) 



Therefore, M{ps{t)pb{t)} satisfies the Liouville equation Eq. (2) and is of course identical 
with the density matrix p{t) of the total system. 

With the help of the Ito calculus, therefore, we are able to illustrate how the interaction 
between the system and bath is decoupled, as shown in Eqs. (3) and (4). As a consequence, 
the system as well as the bath is subjected to random fields. Again, it should be stressed 
that we want to calculate the reduced density matrix, ps{t) — 'TtbM{ps{t)pb{t)}. To this 



end we change the operation order of the tracing over the degrees of freedom of the bath 
and the stochastic averaging to obtain 

Ps{t) = M{psit)Tnpb{t)}. (6) 

The formal expression of TTbPf,{t) can be feasibly acquired by solving its equation of motion 
Eq. (4), 

TrhPbit) = exp I ^ ^ dt' - ip4a{t')] ga{t') I , (7) 

where introduced are the bath-induced stochastic fields 

ga{t) = Tn{gaPbit)} (8) 

with Pb{t) = Pb{t) /T^^'^bPhit) the normalized density matrix for the random bath. Now comes 
a crucial step. We can change the measure of stochastic processes to absorb the trace of the 
bath, that is, we modify the Wiener process Wia{t) asWiait) + 2 jldt'g^{t')/Vh. Complying 
with this, ps{t) will also change accordingly. The mathematical manipulation is nothing but 

571]. To illustrate this transformation clearly, we again 



the Girsanov transformation 



56 



resort to the discrete-time representation. Now inserting Eq. (7) into Eq. (6), one has 

ps{t)= I n n n '^^^oMP {pjm) ^^p < ^ 5Z ^^^^m - ^Pom] Mtj) \ ps it, {pjm}) ■ 

a k=l j [v/i „ J 

(9) 

Note that in this expression one can put the contribution due to the bath-induced field into 
the distribution function P (/ij,i(4)a) , namely, 

n PiPj,la)P{PjAa) ] 7^ [/^i.l" ~ iPjAa] ga{tj) \ 

a I. * Q J 

= (Pj, la - -^(jaitj)^ P (^Pj,4a + -^ga{tj)^ ■ 

Now define new variables pj^ia = Pj^ia — gaitj)/\/h, Pj,2a = Pj,2a, Pj,3a = Pj,3a, and pj^^a = 
Pj,4:a + Waitj)/ \/h. Because both Pj^ia and independent of ga{tj), the Jacobian of 

the variable change is equal to one. With these new variables, therefore, Eq. (9) becomes 

Psit) = j n n n ^f^J^P {P^m) Ps (^t, j/^jMa + -^ga{tj),Pj,2a, Pj,3a, Pj,4a " "^^"'^^J'^}) 
= M{ps [t, {Pka}]}, 
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where the second equahty results from the continuous representation. Notice that ps{t) = 
Ps[t, {/J'ka}] satisfies Eq. (3). Changing the underlying random processes in Eq. (3) 
from {fj^ka} to {flka}, one obtains the equation for the new density matrix ps{t) = 

Ps [t, {ha}], 



For brevity, here and in the following, the functional dependence of the random density 
matrix Ps(t) on {pka} will not indicated explicitly. 

Given ga{t), one only needs to solve stochastic differential equation Eq. (10) and calculate 
the random average to obtain the exact reduced density matrix Ps{t)- The bath-induced 
field in general evoke fast motion and thus ps{t) as well as the white noises to which Ps{t) is 
subjected displays a smaller time scale than its average Ps{t). In the following section we will 
consider the Caldeira-Lcggett (CL) model for which the bath induced field can be worked 
out explicitly. We would like to emphasize that in the CL model there is only one interaction 
term, while within the stochastic description one can deal with many interaction terms. 

B. Caldeira-Leggett Model 

The Hamiltonian of a CL model reads 



where the first two terms define the Hamiltonian of the system and the third one induces 
the Hamiltonian of the bath, the interaction, and a counter-term. The latter can directly 
be read off as 





(11) 




where 




(12) 



J{ijj) is the spectral density function 
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(13) 
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which exactly characterizes the effect of the environment. In general the system-bath linear 
coupling not only provide a dissipation mechanism for the system but also tend to renor- 
malize the potential V{x). To compensate for this renormalization effect the counter-term is 
introduced. It ensures that the system cannot lower its potential energy below the original 
uncoupled value 

Assume that the initial state of the bath is in a thermal equilibrium state for the non- 
interacting harmonic oscillators, that is, pb{0) = e~'^^'' /TTh{e~'^^''}. Because of the linearity, 
the dynamics of the random bath described by Eq. (4) is analytically solvable. Moreover, 
since there are no interactions among the bath modes, the evolution of the bath is fully 
determined by that of the individual mode. In other words, the solution of Eq. (4) can be 
written as 

Pb{t) = llu,,iit,0)phi0)l[u,,2i0,t), (14) 

j j 

where 0) is the forward propagator of the mode j dictated by the Hamiltonian 



and Mj^2(0,t) is the backward propagator dictated by 

1 .2 2^ , 



with the stochastic fields z^±(t) = fJ'2{t) + ipsit) ± ifii(t) ± Pi{t). For a driven harmonic 
oscillator 

h{t) = ^ + 2"^^o^^ + '^(^)^' (15) 
its corresponding propagator u{t, 0) can be worked out in terms of the operator algebra 
method [58y. To do this, it would be better to work in the second quantization formalism. 
Now introduce the creation and annihilation operators. 



2h ^ ^/mmjf ^^^^ 



and 



2h y/2hmujQ 

From the commutation relation [x,p] = ih, one readily finds [a, a^] = 1. The Hamiltonian 
Eq. (15) then becomes 

Mt) = /.o + y^Kt) (a + at), 
8 



where ho = huo (^a^a + l/2j is the Hamiltonian of an isolated harmonic oscillator. The 
propagator un(t,0 ) d ictated by fiQ has been well known and the result can be found in 



several textbooks 
have 



24 



59l-l6l| . For the driven case one uses the interaction representation to 



h 



it) (ae-''^°* + aV"^"*) Ui{t,0) 



dt V 2muJo 

with the initial condition m/(0,0) = 1. Resorting to the Baker-Campbell-HausdorfF for- 
mula 



62j, we readily obtain 

0) = exp [(i(t)a] exp [e(t)a^] exp[j(t)], 



(18) 



where 



and 



d{t) 
e{t) 

m 



y/2mhuio Jo 



V 2mf)wo Jo 



dti exp {—iojoti) i^iti] 



dti exp (ioJoti) z/(ti 



— / dti dt2u{ti)u{t2) exp[iuJo{ti - t2)]. 
2mnwo Jo Jo 



Upon using Eqs. (16) and (17) and Baker-Campbell-Hausdorff formula again, the propaga- 
tor uj{t,0) becomes 



uj{t, 0) = exp 



■mojQ 
~2h 



[d{t) + e(t)] X > exp 



^W)-eW]p}exp{l[A*)-e=W] 



Finally with the given expressions of Uo(t,0) and u/(t, 0) the propagator M(t, 0) = 
Mo(t, 0)M/(t, 0) is also available. Therefore, we propose a straightforward procedure to cal- 
culate Mj,i(t,0) and ■Uj,2(i^,0). Inserting these results into Eq. (8) and calculating the trace 
we finally find 

git) = Vh [ dt' {aR{t - t') [/ii(t') - tfi^it')] + aiit - t') [/i2(t') + ili0)]} , (19) 

where aRj{t) turn out to be the real and imaginary parts of the autocorrelation function of 
the "force" g = Ylj^j^j^ namely, 

1 



a{t) = Ti,{p,{Q)g{t)m} 



duJ{uj) 



coth cos ut — I sm ut 



(20) 



with g{t) = e^Htt/h^^-iHtt/h^ ^^^^^ that aji{t) is dependent on temperature, reflecting 
dissipation effects contributed from both quantum and "classical" motion, while ai(t) is 
independent of temperature, embodying a pure quantum effect. 



C. Formal Solution 



As one can understand, we start from Eq. (10), the working equation in the framework 
of stochastic description. The formal solution of this linear equation is 



Ps{t) = Ui{t,0)ps{0)U2{0,t), 



(21) 



where Ui{t, 0) is the forward propagator of the stochastic system associated with the Hamil- 
tonian 



9a{t) + —V+a{t) 



fa 



(22) 



while U2{0,t) is the backward propagator associated with the Hamiltonian 



E 



fa 



(23) 



with rj±a(t) = Plait) + ipiaif) ± ip2ait) ± Psait) . When we can take stochastic average 
of random density matrix Ps{t) in principle, we can obtain the exact reduced density ma- 
trix Ps{t)- To derive the equation of motion for Ps{t) or the master equation, however, it 
would be better to start from Eq. (10). Whenever the statistical average of the right hand 
side of Eq. (10) can be expressed explicitly in terms of Ps{t) and other known operators of 
the system, one obtains a master equation. 

To acquire the deterministic equation from the corresponding stochastic differential equa- 
tion for arbitrary noises is a challenging, if not impossible task 63l-l66| . In the literature, one 



'or Gaussian 



67| 



may find incorrect results about the derivation of the deterministic equation 
noises and similar statements in developing master equation for open systems 

For the Caldeira-Leggett model, upon substituting the expression of the bath induced 
fields Eq. (19) and carrying out the stochastic averaging, Eq. (10) becomes 



ih 



dt 



Hs,psit) 
t 



+ 



dt'aR{t-t')M{ps,i{t,t')} 



+ 



dt'aj{t-t')M{ps,2{t,t')} 







(24) 



where ps,i{t,t') = \/h[pi{t') - ipiit')] psit) and ps,2{t,t') = \fh[p2{t') + ips{t')] ps{t). In the 
above derivation the nonanticipating property of Ps(t), namely, M {ps{t)dWi(2)it)} = is 
used. It is obvious that when there is no dissipation, only the first term remains on the right 
hand side of Eq. (24). Of course, this is the Liouville equation of the system as it should 
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be. As the stochastic fields are added, the other two terms naturally come out. These terms 
lead to the irreversibility caused by the coupling to the bath d la its induced random fields. 
To obtain the expressions of M{ps_i(2)(t, t')} we resort to the Furutsu-Novikov theorem |68| . 
that is, 

M{Mm^]} = M|^| (25) 

for a white noise n{t) and its arbitrary functional F[fi]. Using this theorem allows us to 
write down 

m{p.,(,0}^^m{||1-,|j|}.o.,(,o 

and 

m{p.,(m')}^v/Sm{|§ + >^}-6.,(m'). 

Therefore, M {p,,i(2)(t, t')} in Eq. (24) can be replaced respectively by Os^i{2){t,t'). For 
brevity, they are called the dissipation operators. With the formal solution of Ps(t) Eq. (21), 
one can directly work out the expression of the functional derivatives 

a,i(t,t') =-^M|a;i(t,tV(t) -p,(t)x2(t,t')} (26) 

and 

a,2(t,t') =M \xiit,t')psit) + psit)x2it,t)} , (27) 



where Xi(t,t') = Ui(t,t')xUi(t' ,t) and X2(t,t') = U2(t,t')xU2(t' ,t) are the Heisenberg op- 
erators. To build up a master equation, of course, we still need to find the expressions 
of Os^i(^2)(t,t') in terms of Ps(t) and other known operators. 

Some comments are in orders. For arbitrary response functions (t), deriving the mas- 
ter equation relies on whether the explicit expressions of dissipative operators Os^i{2){t, t') can 
be worked out. However, when a/j(/)(t) are localized distribution functions, say <yR(i){t) = 
C(R{i)S(t), Eq- (24) becomes 

i^^^^ = [Hs, Ps{t)] + [x, aRM{ps,i{t, t)}] + [x, aiM{ps,2{t, t)}] . (28) 

From Eqs. (26) and (27), we know M{p,,i(t,t)} = -i[x,p,(t)] and M {ps,2it,t)} = 
{x,ps{t)}. Inserting into Eq. (28) yields the master equation 

= Hs,ps{t) -iaji[x,[x,ps{t)]]+ai[x,{x,ps{t)}]. (29) 

Therefore, both the specificities of the system and the bath-induced field determine the 



"existence" of the master equation. 
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III. DERIVATION OF MASTER EQUATION 



Now we shall derive the master equation of a dissipative harmonic oscillator described by 
the Caldeira-Leggett model with V{x) — MuJqx'^/2. As discussed above, the system-bath 
coupling modifies the potential of the system by addition of a counter-term. In the present 
case the frequency of the harmonic oscillator becomes a renormalized one, 

2 2 I ~2 

U — Uq +UJ . 

To acquire the master equation we use the equation of motion in the Heisenberg represen- 
tation to find the operators Xi(^2){t,t') in Eqs. (26) and (27). After some algebra we obtain 

ft 



Xi(2)[t,t ) —cosuj[t — t )x — P ~ / "^isma;(ii— t 

Inserting into Eqs. (26) and (27) and carrying out the involved stochastic averaging, we find 



^(^i) + ^V±{ti) 



Os,i{t,t') = -icosu}{t-t') [x,ps{t)] + sin uj{t-t') \p,ps{t)] 

^ ilL /, "^^^ / "^^^ "^^^^ ~ t')ai(h - t2)ds,i{t, t2) (30) 



and 



A / /X / /X ~ / M SmUj(t — t') ~ / M 

0,,2{t, t') = COSUj{t - t') {X, ps{t)} ' {p, ps{t)} 



2 



Muj 
2 



/ dti I dt2smuj{ti-t')aR{ti-t2)ds,i{t,t2) 
Jt' Jo 

/ dh dt2sinuj{ti-t')ai{ti-t2)ds,2{t,t2). (31) 
Jt' Jo 



-t ftl 

In the derivation we have used the following functional derivatives 



and 

ft 



which are also readily obtained with the formal solution of ps{t) and the Furutsu-Novikov 
theorem. For the integral equations (30) and (31) one can employ iteration to show that 
their solutions Os,i{2){t-,t') possess the following forms 

6s,i{t,t') ^xn{t,t') \p,ps{t)]+xi2{t,t') [x,ps{t)] (32) 
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and 



a,2(t, t') =X2l{t, t) {£, Psit)} + X22(t, t'){p, Psit)} + X2^{t, t') [p, Psif)] 

+ X2A{t,t')[x,Ps{t)]. (33) 

Because the operators [x^psit)]^ [PiPsit)]^ and {p.psit)} are arbitrary, their coef- 

ficients Xjk{t,t') in Eqs. (32) and (33) are determined by Eqs. (30) and (31), which satisfy 
the following integral equations 

xn{t,t') 



-^sina;(t - t') + y dh J dt2smu}{ti - t')ai{ti - t2)xii{t,t2) 
Xi2{t, t') ^ -i cos uj{t - t') + I dti I dt2smoj{ti - t')ai{ti - t2)xi2{t,t2), 

J-ViUJ Jt' Jti 

X2i{t, t') ^ COS uj{t - t') - I dti I dt2smuj{ti-t')ai{ti-t2)x2i{t,t2), 

l-Vii^ Jt' Jo 

X22(t,t') -■ 
X23{t,t') ■ 



sinuit-f) 2 
2 



dti / dt2smuj(ti - t')ai{ti - t2)x22(t,t2), 
Jt' Jo 



(34) 
(35) 
(36) 
(37) 



Muj 



/ dti sin Lu {ti - t') / dt2aR{ti - t2)xii{t,t2) + / dt2ai{ti - t2)x23{t,t2 
Jt' Uo Jo 

(38) 



and 

X24{t,t') 



l-t r l-t nti 

/ dtismu{ti - t') / dt2aR{ti - t2)xi2{t,t2) + I dt2ai{ti - t2)x24{t,t2) 
Jt' Uo Jo 

(39) 



With these expressions at hand, Eq. (24) immediately becomes the desired master equation, 
namely, 

dps{t) 



ih- 



dt 



Hs, p.(t)J + Ai{t) [x, {£, ps{t)}] + A2{t) [x, {p, ps{t)}] 

+ A-i{t) [x, \p, p,{t)]] + Mt) [x, [x, Ps{t)]] , 

where Hs = p^ j (2M) + MuP'x?' jl and the coefficients are 

Aiit) = / dt'ai{t-t')x2iit,t'), 
Jo 

A2{t) = / dt'ai{t-t')x22{t,t'), 
Jo 

^3(t)= / dt'[aR{t-t')xii{t,t') + aj{t-t')x23{t,t')], 
Jo 



(40) 

(41) 
(42) 
(43) 
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and 



Ai{t)= [ dt'[aii{t-t')xu{t,t') + ai{t-t')x2A{t,t')] 
Jo 

As clearly clarified in the literature 



(44) 



38 



Ai{t) is the coefficient for the frequency shift 
because [x, {x, ps(t)}] = [x'^,ps{t)], A2{t) is a quantum dissipation term, A2,{t) reflects the 
anomalous quantum diffusion while A^lt) is the coefficient for the normal quantum diffusion. 



IV. EQUIVALENT TO THE HU-PAZ-ZHANG EQUATION 



The master equation and the dynamics of the dissipative harmonic oscillator has been 
derived and studied by several researchers with diversified methods Dekker used the 

canonical quantization method [35| , while Haake and Reibold employed the Wigner function 
method. The latter also studied the low-temperature and strong-damping anomalies 36|. 
Grabert, Schramm, and Ingold went beyond the factorized initial condition js^]- It would 
be stressed that the master equation in a general environment was derived by Hu, Paz, and 
Zhang (HPZ) by virtue of path integral technique 38|]. An elementary derivation of HPZ 
equation was documented by Halliwell and Yu with the Wigner function approach 39|] . Kar- 
rlein and Grabert pointed out that there is an exact dissipative Liouville operator for certain 
correlated initial conditions 4^. The master equation was also considered by^ Calzetta et 



al. through a stochastic method based on the quantum Langevin equation 41)]. Ford and 



O'Connell 



42| rederived the HPZ equation with the quantum Langevin equation method 



and also analyzed its solution. Strunz and Yu offered an alternative derivation, using the 



quantum trajectory method 



43|. The HPZ master equation has recently been used by Chou, 



Yu, and Hu to derive the master equation of two and more coupled harmonic oscillators in a 



bosonic bath 



44 



451]. Furthermore, the HPZ master equation was extended by Xu et al. 46] 



to the case where time dependent external fields are applied. 

For convenience of comparison, we will always use the results in the paper by Halliwell 
and Yu jsO^. The HPZ master equation takes the same form as Eq. (40). The corresponding 
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coefficients in our notation read 



and 



wliere 



(45) 



Bi{t) = f dt'ai{t-t')x2i{t,t'), 
Jo 

B2{t)= I dt'ai{t-t')x22{t,t'), (46) 
Jo 

Bsit)^ I dt'[aR{t-t')xu{t,t') + ai{t-t')x23{t,t')], (47) 
Jo 



,{t) = [ dt' [aR{t - t')xi2{t, t') + ai{t - t')x24{t, t')] , (48) 
Jo 



X2l{t,t')=U2{t')-^U,{t'), (49) 

xi,{t,t')^^G,{t,t'), (51) 

xi2{t,t') = -tG[{t,t'), (52) 

X23{t,t') ^ — J dti J dt2aR{ti-t2)Gi{t,t2)G2{t',h), (53) 



and 



X24it, t') = c^ti ^ di2aH(*i - ia)^^*, i2)G'2(t', ^i). (54) 

Here tfie dot over Uj{t) {j — 1, 2) stands for tfie derivative witfi respect to t and tfie func- 
tions Uj{t) are tfie soiutions of tfie fiomogeneous integro-differentiai equation 

+ <t) + A £ dt'aiit - t')u{t') = (55) 

with inhomogeneous boundary conditions 'Ui(O) = 1, Ui(t) = and ^2(0) = 0, U2{t) = 
1 and Gj{ti,t2) {j = 1,2) are tfie Green's functions obeying 



V dif ^ J ^^^^MJ "^^^"^^^1 ~ ^3)G(t3, ^2) = 5(tl - t2) 



(56) 



witfi specified initiai conditions at tfie fixed initiaf and finaf times Gi{ti — 0, ^2) = 0, G[{ti — 
0,^2) = and 6^2(^1 = ^,^2) = 0, G2{ti — ^2) = 0. Here tfie prime in G'j{ti,t2) stands for 
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the derivative with respect to the first variable, that is, G'j{ti, ^2) = dGj{ti, t2)/dti. Because 
of causahty Gi{ti, ^2) = for ti < ^2, while (^2(^1, ^2) = for ti > ^2- 

Now we show that the HPZ equation and that derived with the stochastic formulation 
are identical. To this end we only need to prove that Aj{t) = Bj{t) {j = 1 — 4), respec- 
tively. As clearly shown in Eqs. (41) — (44) and Eqs. (45) — (48), all functions Aj{t) and the 
counterparts Bj{t) are integrals over the time range [0,t]. Therefore, a sufficient condition 
for Aj{t) — Bj{t) is that the corresponding integrands are the same. Besides, because these 
integrands consist of factors aR(^i) (t) that are dependent on the specificity of the dissipation 
and can be arbitrary, one can further simplify the proof significantly. 

A. Proof of Ai (t) = Bi (t) , A2 (t) = B2 (t) 

To prove Ai{t) = Bi{t), one should prove X2i{t,t') = X2i{t,t'). Note that Uj{t) satisfy 
linear differential Eq. (55). Because X2i{t,t') is a linear combination of ui{t') and U2{t'), as 
a function of t', it should also obey Eq. (55), 

i ^ + ^2iit, t') ^-j^ J dhaiit' - h)x2iit, h). (57) 

On the other hand, from the integral Eq. (36) we can show by a straightforward algebra that 
calculating the second order derivative of X2i{t, t') on both sides one can obtain the equation 
the same as Eq. (57). Moreover, X2i{t,t') \t'=t= X2i{t,t') \t'=t— 1 and dx2i{t,t')/dt' \t'=t— 
dx2i{t , t') / dt' \t'=t= 0, the initial conditions are the same. Therefore, Ai{t) is identical 
with Bi{t). On the same line, one can prove A2{t) = B2{t). 

B. Proof of A3{t) = Bz{t), Ai{t) = Bi{t) 

As pointed out above, the problem of proving A2,{t) — B2,{t) can be changed to proving 
the equivalence of the two involved integrands. That is, one needs to show xii{t,t') — 

xn{t,t') and X2z{t,t') = X23{t,t'). 

Let us look at Eq. (56) with the initial condition 6*1(0, 12) = and G[{0, 12) = 0. Suppose 
the second term on the left hand side is given. Then Eq. (56) can be viewed as a function 
of ti and can be "solved" with Green's function method. Now the required Green's function 
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/ dt3 / dt4SmUj{ti - t3)ai{t3 - t4)Gi{t4,t2). 

Jt2 Jt2 



obeys 

{§1+^') G{t„T)^S{t,-T) 

with G{ti,T) \ti<T= and dG{ti,T)/dti |ti<T= 0. Its solution is 

^, , sina;(ii — r) ^ , 
G{h,T) = ^— '-9{ti - r). 

Therefore, the "solution" of Eq. (56) can be written as 

Giih, t2) = J drG{h,T) 5(t -h)-^ j dUai{T - U)G^{U, ^2) 

_sinw(ti-t2) 2 '•'^ 
~ uj Mcu 

Substituting into Eq. (51) yields 

xii(t,t') = -^sini^(t -t') - dh J dt2smu{t-ti)ai{ti-t2)xii{t2,t'). (58) 

Note that from Eq. (20) one can see Q;/(ii — ^2) = —oii{t2—ti). Making change of integration 
orders and variables in the double integral then leads to 

^ii(^,i') = j^^sina;(t - i') + dh J dt2smuj{t - t2)ai{ti - t2)xii{ti,t'), (59) 

which is an integral equation of Xii{t, t') with respect to the first argument t. Now we show 
that Xu{t,t') determined by Eq. (59) is identical with xu{t,t') solved from Eq. (34) that 
is an integral equation with respect to the second argument t'. To this end, we discretize 
the variables t and t' so that Xii{t,t') and xii{t,t') can be represented as matrices. To be 
specific, the elements of matrices Xq, X, X, and a take the discretized values of xo{t,t'), 
xii{t,t'), xii(t,t'), and a(t — t'), respectively. Here xo(t,t') = sina;(t — t') is introduced. 
With these matrices, the integrals in Eqs. (34) and (59) become matrix products, 

X^Xo-XaXo (60) 

and 

JT = Xo - XqCxX. (61) 
Solving these matrix equation with elementary algebraic manipulations, we find 

X^Xo{l + cxXo)-^ (62) 
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and 



X' = {1 + Xqcx)-^ Xo. (63) 

Because (1 + XqCx.) Xq = Xq (1 + cx.Xq), one immediately obtains X = X'. Therefore, 
we find Xu{t,t') = xu(t,t'). Thus, Eq. (51) can be recast as Giit,t') = —iMxnit.t'). 
Substituting into Eq. (53) leads to 



X23{t,t') = ^ j^^ Jo ^^2"^*^^^ -h)xii{t,t2)G2{t',ti] 



(64) 



We now treat (^2(^1, ^2) in the same way as we did for 6*1(^1, t2) in the above. As a result, 
^2(^1,^2) satisfies the following integral equation, 

G2{ti,t2) = ^^^^('^^-'^2) _ J_ /" r dUsinuj{h-t{)ai{h-U)G2{UM)- (65) 

Inserting into Eq. (64) and rearranging, we find the integral equation of X23{t,t') is iden- 
tical with that of X23{t,t'), Eq. (38). Therefore, X23(t,t') = X23(t,t'). Similarly, we can 
demonstrate the equality A4{t) = Bi{t). 



V. MASTER EQUATION OF DRIVEN HARMONIC OSCILLATOR 

Consider a dissipative harmonic oscillator driven by general external time-dependent 



fields 46144 9(1 with Hamiltonian 



Hs{t) = ^ + \mu:'x' + h{t)x + /2(t)p. 

We shall work out the master equation in the driving case along the same line of deriving the 
master equation of the dissipative harmonic oscillator in Sec. III. Note that for the external 
time- dependent fields only act on the system, the bath-induced field is the same as that 
without driving fields. We can thus start from Eqs. (26) and (27). By solving the equations 
of motion for the Heisenberg operators and taking stochastic averaging, we find 

Os,i{t,t') = -icosuj{t-t') [x,ps{t)] + — —sin u{t -t') [p,ps{t)] 

IVI LO 

-t rt 



^Wuj It, '^^^ It ^^2sinw(ti - t')a/(^i - ^2)0^,1(^,^2) 
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and 



ds,2{t,t') ^COSUj{t - t') {X,ps{t)} - ^'"""^^^ {p,Ps{t)} 



2 



Muj .h, 
2 



Muj Jh 



Muj 



/ dh dt2smuj{ti - t')aR{ti - t2)ds,i{t,t2) 
Jt' Jo 

/ dti / dt2smuj{ti - t')ai{ti - t2)ds,2it,t2) 
Jt' Jo 

[ dtisinuj{ti-t')fi{h)ps{t)-2[ dtiC0SUj{h-t')f2{ti)ps{t). 
Jt' Jt' 



't' Jt' 
Similar to the undriven case, these dissipation operators possess the following forms 

ds,i{t, t') =xn{t, t') IP, ps{t)] + xu{t, t') [x, ps{t)] 
and 

Os,2{t, t') =X2l{t, t') {X, psit)} + X22{t, t') {p, Psit)} + X2z{t, t') [p, + X2A{t, t') [£, 

+ X2->{t,t')ps{t), 

where all coefficients Xjk{t,t') except X2^{t,t') are the same as that of undriven case 
Eqs. (34) — (39). The new term is determined by 

a;25(t,t') = - J^^^ dtismuj{ti-t')fi{ti)-2 J^^ dh cos uj{ti - t')f2{ti) 

— / "^^1 / dt2SYilUj{li,\-i!)ai(^i-t2)X2hi^-,t2). 

Muj J^, Jq 



Therefore, substituting M{ps^i(2){t, t')} by Os,i(2){t, t') in Eq. (24) yields the required master 
equation 

^^dpM 



dt 



Heffit),psit)\ +Ai{t) [X,{x,psit)}]+A2it) [X,{p,ps{t)}] 

+ A^it) [x, [p, ps{t)]] + A^{t) [x, [x, ps{t)]] , 
where the effective Hamiltonian H^ffit) is 

4//(*) = Hs{t) +x [ dt'aiit - t')x25it, t'). 
Jo 

Here again, the coefficients Aj{t) {j — 1 — 4) are the same as the undriven case, which are 
given through Eqs. (41)— (44). It is obvious that the second term in Heff{t) reflects the very 
interplay between the system and the bath mediated by the driving fields. 
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VI. CONCLUSION 



Classical dynamics of dissipative systems is traditionally described by the Langevin equa- 
tion. It has been shown that the stochastic formulation provides a similar description to 



quantum dissipative systems [28|. In this framework the dissipative system obeys Liou- 
ville equations subjected to the stochastic fields due to the bath. Based on the stochastic 
formulation, flexible numerical methods have been proposed and used to solve dissipative dy- 
namics 



30l |. As complementing to previous work the present paper provides conceptual and 
analytical results. We have illustrated how to obtain the bath-induced field of the Caldeira- 
Leggett model through elementary solution of quantum linear systems. Furthermore, we 
have elaborated a systematic approach to derive the master equation, if it exists. 

The reduced density matrix can in principle be obtained by solving the stochastic Liou- 
ville equation and calculating stochastic average. Because it is difficult to have convergent 
stochastic averaging, a master equation describing the evolution of the reduced density ma- 
trix is highly desired. We have shown the existence of the master equation relies not only 
on the feature of the dissipation charactered by the spectral density function, but also on 
the dynamics of the stochastic system itself. For linear systems, we have found that the 
"dissipative operator" due to the interplay of the system and the stochastic field is exactly 
solvable and thereby derived the master equation. We have shown that the master equa- 
tion is equivalent to the HPZ equation derived by Hu, Paz, and Zhang using path-integral 



approach 



38, 



39[. We would like to point out that the coefficients in our master equation 



are determined by a set of integral equations, which may not suffer from the mathematical 



50| 



problems concerned by Fleming, Roura, and Hu 

We have also shown that the master equation of driven harmonic oscillator can be derived 
similarly. In this case, the system is dressed by both the driving and the stochastic fields, 
although the dissipation operators appear the same as that of the undriven case. 
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